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Abstract 

Several issues related to pentaquark searches relevant for current and fu- 
ture high-energy colliding experiments are discussed. We make an attempt to 
explain why pentaquark candidates are not seen by some experiments, and 
what makes the HERA experiments so special in such searches. 



1 Introduction 



Recently, many experimental groups made significant efforts to find new states which 
can be explained as consisting of five quarks. A number of experiments [1,2] in- 
cluding ZEUS [3] have reported narrow signals in the vicinity of 1530 MeV in the 
nK + and K s p invariant-mass spectra. The signals were consistent with the exotic 
pentaquark baryon state + (153O) with quark content uudds [4]. Several other ex- 
periments searched for this state with negative results [5-9]. ZEUS also reported 
negative searches for other pentaquark candidates [10-12]. 

At present, ZEUS is the only high-energy colliding experiment which observes 
the G + (1530) state in the Kg p decay channel. The HI experiment has evidence 
for the lightest charmed pentaquark state [13]. Several other colliding experiments, 
BaBar [8] , ALEPH [6] , Belle [9] and CDF [7] do not observe any of the pentaquarks. 
Although in some cases the results are still preliminary, it is evident that the statis- 
tics of such experiments for several known states exceed the statistics of those ex- 
periments which observe the pentaquarks. While non-observation of pentaquarks 
in e + e" collisions does not necessarily contradict to the observation of a signal pre- 
dominantly produced in the forward region with a baryon in the beam or in the 
target, it is more difficult to explain why no pentaquarks are seen in pp collisions at 
Tevatron and for some other fixed-target experiments. 

The aim of this article is not to give a detailed review of the present status 
of pentaquark searches, which can be found elsewhere [14,15], but to discuss the 
most intriguing problem: why the ZEUS experiment observes the O + (1530) state 
in deep inelastic scattering (DIS) at medium Q 2 . The answer on this question may 
also explain the negative searches at Tevatron, and perhaps, the negative results 
for some fixed-target experiments. We also will discuss several experimental issues 
which could be important for current and future pentaquark searches in high-energy 
experiments. 

The existing experimental results are still rather controversial, therefore, it is 
difficult to discuss them without making any assumptions. We assume that the 
+ (153O) state in the K s p decay mode is indeed a pentaquark, but not a new 
£ state with a strongly suppressed An decay mode. Secondly, at present, there 
is rather weak evidence for antipentaquarks: the ZEUS result [3] indicates that 
the peak at 1522 MeV is mainly driven by the K s p decay mode rather than due 
to the K s p channel (although a structure at 1480 MeV, which is consistent with 
the PDG £(1480) bump, is more pronounced for the K s p invariant mass). Since 
antipentaquarks can only be produced by soft fragmentation, the second observation 
avoids the contradiction with negative searches in e + e~ annihilation experiments. 

2 Possible production mechanism 

We will start with the discussion of possible differences in the production mechanism 
of baryonic states in ep collisions at HERA and other experiments. 
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Probably, the most simple explanation of why pentaquarks have not been ob- 
served in e + e~ experiments lays in the fact that these experiments do not have the 
net baryon number in the colliding beams. If one thinks in terms of the pentaquark 
explanation of the O + (1530) signal, the production of this state in e + e~ should be 
significantly suppressed since it requires five quarks to be glued together by the soft 
fragmentation mechanism - something which is difficult to explain using the conven- 
tional hadronisation mechanisms. This could explain the negative results in e + e~ 
annihilation experiments. 

The situation is different when a baryonic state is present in either beam and/or 
in the target. In this case, the fragmentation mechanism is not as simple as in 
e + e~, and it is possible that the hadronic final state would have enhanced baryon 
production rate in comparison with the e + e _ annihilation. If the presence of the 
net baryon number is indeed a necessary condition for formation of the + (153O) 
pentaquark, then one should explain the negative results of the CDF experiment [7] . 

It should be emphasised that the 6 + (1530) state was observed in DIS at the 
exchanged photon virtuality Q 2 > 20 GeV 2 [3], while no signal was reported in 
photoproduction (Q 2 ~ 0). The photoproduction events are dominated by the 
boson-gluon fusion (BGF) mechanism, 7*g — > qq, which makes the ep collisions to 
be very similar to pp. Thus the negative result from Tevatron is not in contradiction 
with non-observation of the 6 + (1530) state in photoproduction at HERA. 

Note that the photoproduction data selected by ZEUS are trigger-biased mainly 
by the detection of high jets. Therefore, such data are enriched with the BGF- 
type of events. In contrast, DIS data are free from such bias as they were collected 
by triggering the scattered electron, with a minimum bias on the hadronic final 
state. 

What makes the HERA data so special with respect to exotic baryon searches? 
In ep collisions, a significant fraction of events contain a leading baryon carrying a 
large fraction of the incoming proton momentum [16-19]. There are several models 
describing the leading baryon production. One type of models is based on the QCD 
building blocks, quarks and gluons, and explains such baryons as the result of the 
hadronisation of the spectators from the incoming proton. In the second type of 
models, the leading baryons can be explained by the exchange of virtual particles. 
Finally, such baryons could originate from decays of low-mass proton resonances in 
diffractive interactions. 

Recently, both HI [17] and ZEUS [18] reported a rise in the leading-baryon 
yield with increase of Q 2 from the photoproduction regime to DIS. This could be 
explained by absorptive effects. If one assumes that a simitar effect could also 
contribute to the exotic-baryon production, this may explain why ZEUS does not 
see the 6 + (1530) signal for photoproduction and low Q 2 DIS. This could be a good 
guess, if one proves that the most central fragmentation region, corresponding to 
the pseudorapidity range | rj |< 1.5 used for ZEUS pentaquark searches, can also be 
populated with baryons after fragmentation of the nucleon remnant (but, of course, 
with a significantly smaller probability than that for the most forward 77 region 
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accessed by the leading-baryon calorimeters installed by ZEUS and HI). 

To verify this assumption, Monte Carlo (MC) studies were performed using 
the following models: PYTHIA6.2 [20], HERWIG 6.5 [21], ARIADNE 4.12 [22] 
and LEPTO 6.5 [23]. The physics calculations were performed with the module 
"baryons_DIS.rmc" of the RunMC package [24]. The MC models do not contain 
the simulation of the leading-baryon effects due to the exchange of virtual particles, 
but rather they attempt to describe the effect assuming that the nucleon remnant 
is a diquark composed of the valence quarks as spectators. While PYTHIA, ARI- 
ADNE and LEPTO use the LUND model to describe the hadronisation by stretch- 
ing a colour triplet string between the struck quark and the diquark 1 , the HERWIG 
model uses the cluster model to hadronise such diquarks. In case of LEPTO, the 
soft-colour interaction model is used to simulate the rapidity-gaps events. 

At this point, we are not interested in how realistic such models are in the 
description of the production rate of the leading baryons which carry a large fraction 
of the incoming proton energy [17, 18]. The main question is whether they can 
predict events which have an enhanced production rate of the conventional baryons, 
compared to that of antibaryons for the central fragmentation region where the 
O + (1530) state was seen by ZEUS. 

Figure ^ shows the hadron-antihadron asymmetry separately for mesons and 
baryons as functions of px and rj in the laboratory frame. The pseudorapidity was 
defined as r] = — In (tan |) , where the polar angle, 9, is measured with respect to the 
proton beam direction (i.e. positive rj corresponds to the proton beam direction). 
The Monte Carlo simulations were performed with ARIADNE for Q 2 > 10 GeV 2 
and Q 2 > 1000 GeV 2 . The simulations were done for the proton (E p = 920 GeV) 
and the electron (E e = 27.5 GeV) colliding beams. The asymmetry was defined as 
r(h) = N(h)/N(h), where N(h) (N(h)) is the number of counted hadrons (anti- 
hadrons) and h — m,b, where m and b denote mesons and baryons, respectively. 
The minimum transverse momentum for hadrons was set to 150 Me V as for the 
experimental searches [3]. 

It is interesting to observe that the ratios shown in Fig. ^ are above unity, 
i.e. they indicate the asymmetry. This is not a surprise since the hadron-final 
state is affected by the presence of the valence quarks from the incoming proton 
with the total electric charge +1. A similar effect is also expected for all types of 
hadrons (mesons and baryons) in the current jet [25], and it was found to be stronger 
with increase of Q 2 . The most interesting feature of Fig. ^is that the asymmetry 
is larger for baryons than for mesons, and this tendency increases at low pt and 
positive rj. This could be an indication of the presence of the proton debris after 
the fragmentation process. 

Thus, occasionally, baryons can be kicked out of a proton by the fragmentation 
process to the central pseudorapidity region. The probability of such effect is small, 
but not negligible, especially at high Q 2 . As was discussed above, there are other 

1 For BGF, this picture is modified as tree quarks remain in colour octet after the gluon emission. 
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mechanisms responsible for the leading-baryon production, and they may further 
contribute to the baryon-anibaryon asymmetry in the pseudorapidity region shown 
in Fig. 

To illustrate the presence of the net baryon number from the incoming beam 
better, one can calculate the double ratio, r(b)/r(h). This double ratio has the 
advantage that the trivial effects related to leading hadrons in the current jet are 
removed, thus r(b)/r(h) is more sensitive to the presence of the baryonic enhance- 
ment due to the proton-remnant fragmentation. Fig. |2 illustrates such double ratios 
for ARIADNE, using the points shown in Fig. ^ The result shows a significant 
deviation from unity. The effect is at the level of 2% for 77 ~ 2 at Q 2 > 10 GeV 2 . 
For Q 2 > 1000 GeV 2 , the double ratio can be as big as 7% for the same r] region. 
The observed effect is mainly due to baryons produced at low p T . 

The same studies were also performed with LEPTO and HERWIG (not shown). 
For Figs. □ and H the LEPTO model generated for Q 2 > 10 GeV 2 with the soft 
colour interactions is somewhat above the ARIADNE model for Q 2 > 10 GeV 2 , 
but it is below the ARIADNE predictions for Q 2 > 1000 GeV 2 . HERWIG shows 
a stronger asymmetry than LEPTO for both the ratio and the double ratio, which 
is likely to be related to differences between the Lund and cluster hadronisation 
models, fragmentation mechanism. Thus the ratios shown in Figs. ^ and 121 have 
large mo del- dependent uncertainties. 

Figures El and |U show similar ratios for photoproduction events obtained using 
PYTHIA and HERWIG. Both resolved and direct contributions were included in 
the simulation. The events were generated by requiring the minimum transverse 
momenta E™ m = 8 GeV for the hard subprocess. The asymmetry shown in Fig El is 
significantly smaller than for DIS. For the double ratio, the effect is below 0.5% for 
-2 < 77 < 2. 

Thus, the MC simulations indicate that, for DIS events, a contribution from the 
fragmentation of the proton remnant is present even in the central fragmentation 
region, | 77 |< 1.5, i.e. where the O + (1530) state was observed by ZEUS. This 
observation has a potential to explain the presence of the O + (1530) state as well as 
other possible exotic baryonic states in ZEUS DIS data [3,11]. At the same time, 
this mechanism leads to a reduction of the O + (1530) production rate at low Q 2 and 
photoproduction, which is consistent with the ZEUS studies of O + (1530) peak 2 . 
Since there is no strong evidence for antipentaquarks, the proposed explanation is 
not in contradiction with the ZEUS data. 

The baryon-antibaryon asymmetry can also be studied by reconstructing A/A, 
p/p etc. ratios. For the A/A ratio [11], there is indeed a large asymmetry, but it is 
difficult to decouple it if from secondary-scattering interactions. Recent study of 5 
baryons [11] shows that the current statistical precision of HERA-I data is not yet 
sufficient to establish the H/H asymmetry at the 1 — 2% level. One should again 
note that the baryon-antibaryon ratio measured in the central pseudorapidity region 

2 Another possible explanation of the reduction of + (153O) signal at low Q 2 was discussed 
in [3]. 
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cannot give ultimate proof of the enhanced baryon production due to fragmentation 
of the proton remnant, because a similar asymmetry should also be expected for 
all hadrons (including mesons) due to leading hadrons in the current jet. As was 
discussed before, the double ratio is more appropriate tool to observe the discussed 
effect, since the contribution from leading hadrons in the current jet is removed. 

In summary, if the pentaquark production is related to the proton-remnant frag- 
mentation, then pentaquarks should be best seen at low p? i n the laboratory frame, 
and for positive values of rj. The production of antipentaquarks should be sup- 
pressed relative to pentaquarks. In addition, high-Q 2 region is more favourable, as 
the incoming proton receives a stronger kick from the virtual photon; in the LUND 
fragmentation picture, a colour string stretched between the diquark and the struck 
quark has a higher chance for dragging the diquark system to low values of pseudo- 
rapidity than for low Q 2 events. To verify this mechanism, high statistics data from 
HERA-II are necessary in order to perform more differential studies of the + (153O) 
signal. 

3 Experimental reconstruction 

It is noteworthy that, from the experimental point of view, the decay channel 
+ (153O) — > K s p is not simple since the dE/dx method is often used to identify the 
protons, while Kg are reconstructed using the secondary-vertex algorithms. In this 
approach, the identification of protons is only possible at low-momenta (p < 1.5 GeV 
in ZEUS case). However, in this case, the phase space available for 6 + (1530) decay 
is dramatically reduced as the decay kinematics requires that the momenta of the 
proton candidates should be typically higher than the momenta of K s [26]. Note, 
in colliding experiments, a significant fraction of K s mesons have momenta above 
p ~ 1.5 GeV, which is the ZEUS cut used for the proton identification. 

The HERA data [2,3] indicates that the K s p invariant- mass spectrum has non- 
negligible contributions from not well established PDG E states. Such baryons are 
not implemented in MC simulations, therefore, the MC predictions (scaled to the 
luminosity of the data) are significantly below the data for the K s p spectrum. Such 
sensitivity to the E states may also be used by current and future experiments to 
test the sensitivity to 6 + (1530). 

The resolution for the K s p invariant mass is another central problem for pen- 
taquark searches as the width of the O + (1530) baryon is predicted to be narrow. 
ZEUS is among the experiments which have excellent resolution for this study, 
~ 2.5 MeV [3] (using the PDG mass for K s to improve the resolution for the K s p 
mass distribution). Keeping in mind the difficult background near the 1480 MeV 
mass region due to the PDG E(1480) bump, a structure near the 1520 MeV region 
will be difficult to observe if the resolution would be 3 — 4 times larger. 

As was discussed before, if the contribution of the proton remnant to the cen- 
tral fragmentation region is expected at high Q 2 , then one may expect that the 
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pentaquark production is also enhanced at high Q 2 . However, from the experimen- 
tal point of view, pentaquark studies at very high Q 2 could be affected by high 
combinatorial background. This will be discussed in more details below. 

3.1 Combinatorial background and experimental sensitivity 

At high energies, the combinatorial background for the invariant- mass reconstruction 
scales approximately as < n > 2 , where n is the number of produced tracks. Going to 
a high energy (y/s, jet E T , Q 2 , W, etc.) available for hadroproduction, combinatorics 
for fully-inclusive events is rising due to a large production rate of non-heavy flavour 
hadrons produced by soft fragmentation and resonance decays. For HERA, an 
increase in the track multiplicity and combinatorial background is expected for high- 
er photoproduction events, or for very high Q 2 events in DIS. This should mask 
the peaks from hadrons which may carry the quantum numbers of colliding beams, 
or states which can only be produced after the fragmentation of partons created by 
hard subpro cesses. 

The most relevant example illustrating this point is heavy-flavour mesons or 
baryons. The production of such states by pure fragmentation is significantly sup- 
pressed, but they can be produced in the hard BGF subprocess 7*g — > cc at the 
fragmentation stage. Thus, for high-energy processes (or for high-.&r jets), there is 
a stronger combinatorial background than at low energy, since a significant fraction 
of strange kaons and baryons produced by hadronisation and resonance decays can 
contribute to the relevant invariant-masses distributions. 

This point has been checked by using the same MC samples as those shown in 
Figs.^EJ First, the decay channel D° — > K~n + was reconstructed by combining the 
momenta of two final-state charged particles. It was found that, for the ARIADNE 
event sample generated at Q 2 > lOGeV 2 , the ratio of the signal events over the 
background was about 2.3, while for the photoproduction sample of PYTHIA with 
Ejnin _ g 1 this ratio was 1.5. Such difference is mainly related to a high particle 
multiplicity in photoproduction data than in DIS. In contrast, when 0(1020) mesons 
decaying to K + K~ are reconstructed, the signal-over-background ratios were 1.3 and 

1.2 for DIS and photoproduction, respectively. Thus, the signal-over-background 
ratio for the states which can be produced by pure fragmentation (like 0-mesons) is 
about the same for DIS and photoproduction, while it is significantly different for 
hadrons which have a production mechanism driven by the hard subprocess or by 
other mechanisms unrelated to pure fragmentation and resonance decays. 

In this respect, the statement [14] that experiments with superior statistics for 
known hadrons (such as A(1520)) have better chances to find + (153O) is not fully 
correct for high-energy experiments. The production mechanism of such states is 
soft fragmentation, thus the signal-over-background ratio does not change signifi- 
cantly as a function of the energy available for the hadronic final state. In contrast, 
the production mechanism of the G + (1530) state is unlikely related to pure fragmen- 
tation mechanism, thus the combinatorial background for the + (153O) reconstruc- 
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tion is expected to increase with energy. Therefore, high statistics for reconstructed 
non-heavy flavour mesons and baryons can only indicate the quality of particle 
identification. This is an important issue, but not the most important compared to 
possibly exotic production mechanism of the + (153O) state, and the experimental 
difficulties related to the increase of combinatorial background with energy. 

Note that for the A(1520) baryons, which have a similar decay mode K~p as 
the + (153O) state, in some cases it is also difficult to say about the experimental 
sensitivity to the + (153O), even if one assumes the same production mechanism 
as for the pentaquark state. Indeed, the ZEUS preliminary results [12] showing 
the A(1520) signal at Q 2 > 20 GeV 2 does not look impressive compared to the 
corresponding peak from BaBar [8]. This is due to the fact that the reconstruction 
of A(1520) at ZEUS was done in a significantly smaller allowed phase space than 
that for the + (153O), since K ± mesons were reconstructed by using the dE/dx, 
i.e. only K ± mesons with relatively low momenta, p < 0.8 GeV, were used. In 
contrast, the K$ sample used for the + (153O) searches did not have such limitation, 
since no the dE/dx requirement was used (there is only 20% of Kg in the range 
p < 0.8 GeV). The protons used for the A(1520) reconstruction were also selected 
by using a significantly stronger dE/dx cut. 

It seems that the best known state which can be used to check the experimental 
sensitivity to + (153O), and to compare it with other experiments, is A c (using 
the K s p decay channel). For ep collisions, the A c baryons are produced through 
the BGF mechanism, 7*g — > cc. Thus the cross section for A c is decreasing with 
increase of Bjorken x and Q 2 . Unfortunately, if the G + (1530) production mechanism 
proposed in Sect. 121 is correct, the yield of A c baryons decreases for the kinematic 
regions where the 6 + (1530) cross section is expected to be high. 

3.2 Particle misidentification 

Particle misidentifications in hadron spectroscopy can roughly be divided into the 
following categories: 1) misinterpretation of a peak which may correspond to a 
known state when a different mass assumption for the decay products is used, and 
2) misreconstruction of a state after using the same tracks twice (track splitting). 
We will not discuss here the first problem as this is always a first check to be 
made to exclude possible reflections (which usually never produce a very sharp peak 
consistent with the detector resolution). 

Here we will concentrate on the second problem. It was already pointed out by 
several authors [27] (see other details in [14]) that, for example, using ghost tracks 
associated to A-baryon decays can lead to peaks at 1540 MeV. 

Generally, the experiments which have reported the + (153O) state are well 
aware of the danger related to the ghost tracks. The issue is probably more im- 
portant for high-energy experiments with large combinatorial background under the 
observed peak as it is impossible to scan each event by eye in order to verify which 
track belongs to which decaying state. 
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Note that all experiments which observe the K s p peak claim the reconstructed 
G + (1530) masses somewhere between 1520 and 1530 MeV, so the spurious peak 
shown in [14] is by 20 — 10 MeV shifted from the observed peak. At ZEUS, every- 
thing was done to avoid such double counting by removing A hyperons from the 
consideration and by reducing the background under the K s peak as much as possi- 
ble. The contribution from the A decays was further reduced using the tracks fitted 
to the primary vertex for the proton candidates. 

Finally, ZEUS observes the @ + (1530) peak only for a specific type of events, 
mainly for DIS. It is hard to imagine why the spurious peak appears in DIS only, but 
not for high-i?T photoproduction events which have a high average track multiplicity 
and a larger production rate of A baryons compared to DIS. 

The last observation is probably a good hint for future pentaquark searches: if a 
signal disappears for certain class of events (where it can potentially disappear due 
to some physical reasons, and assuming that such class of events was not preselected 
using track- related quantities), this would likely mean that the observed peak is not 
due to track-related misreconstructions. 

3.3 Conclusions 

The HERA collider is unique and plays a significant role in pentaquark searches. 
This uniqueness is in a possible access to the net baryon number from the collid- 
ing proton beam, even when the most central fragmentation region is used for the 
baryon reconstruction. This conclusion is mainly relevant for DIS, while the effect 
is significantly smaller for gluon-driven photoproduction. If the assumptions made 
in this paper are correct, the pentaquark signal should be best seen in the forward 
region of pseudorapidity, at low p T and at medium or high Q 2 regions. The signal 
from antipentaquarks should be suppressed compared to pentaquarks. 

Another feature of ep collisions at HERA is relatively low energy for the hadronic 
final state compared to the Tevatron experiments. This has a clear advantage as it 
leads to a low combinatorial background for the K s p mass reconstruction. Further- 
more, according to several predictions [28], the 6 + (1530) production cross section 
can be suppressed at high energies. 

For the nearest future, HERA could be the only high-energy colliding experiment 
which is well suited for the exotic baryon searches, since the pentaquark production 
in e + e~ annihilation is expected to be suppressed compared to the experiments with 
baryons in the colliding beams. Therefore, the future Linear Collider is likely to be 
not well suited for such searches. From the other hand, the hadronic final state in pp 
collisions at LHC is expected to be too complicated, and a significant effort should 
be made to deal with the combinatorial background for the Kg p invariant mass. 
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Figure 1: The ratio of the total number of final-state hadrons to antihadrons, 
N(h)/N(h), as functions of p T and rj in DIS. The ARIADNE Monte Carlo model 
is used. The ratios are shown for mesons and baryons separately for two regions in 
Q 2 - 



11 



I, 

z 

? 

z 



IE 



\n 







I 


I I 


I I 


I 


I - 


1.04 


















• ARIADNE Q 2 > 


10 GeV 2 




z 


1 .Uo 












~z 


1.02 




O ARIADNE Q 2 > 


1000 GeV 2 




z 












— 


1.01 


















9 








I z 


I 






• 


• 




• — 

I : 










• 


• 




0.99 














0.98 














0.97 








o 




-z\ 


0.96 








o 
















o 


O : 


0.95 




i i i i i 


i i I i i i i I 


I I I I I I I I I I I 


I I I I 


i i i I i i i i _ 






0.5 


1 1.5 


2 2.5 


3 


3.5 A 














P T (GeV) 


1.1 
















— r- 


i i i , i 


i i I i i i i , 


i i i i , i i i i I i i 


, , I , 


, , , i , , , ,_ 


1.08 






















O 


o 


- 

o 


1.06 














1.04 














1.02 







O 


o 


• 


• _: 


1 




• 


• 


• • 






0.98 














0.96 


















I , I I , 


i i I i , i i I 


, I I , I I I , I I I I 


I I I I 





-1.5 



-0.5 



0.5 



1.5 



Figure 2: The double ratio (see the text) as functions of pt and rj for ARIADNE. 
The ratios are shown for two regions in Q 2 . 
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Figure 3: The ratio of the total number of final-state hadrons to antihadrons as 
functions of pt and rj. The photoproduction events were generated with PYTHIA 
6.2 and HERWIG 6.5 using the minimum transverse momentum for the hard scale 
E™ in = 8 GeV. The ratios are shown for mesons and baryons separately. 
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Figure 4: The double ratio (see the text) as functions of px and i] for photoproduc- 
tion using the minimum transverse momentum for the hard scale E™ m = 8 GeV. 
Monte Carlo studies were performed using PYTHIA 6.2 and HERWIG 6.5. 
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